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SECTION  I 
SUMMARY 


The  purpose  of  this  phase  of  the  Ballistic  Research  Labora¬ 
tory  contract  was  investigation  of  the  feasibility  of  developing 
a  simple  analytical  model  which  could  be  used  to  predict  loads  on 
a  penetrator  in  semi-infinite  concrete.  In  its  ultimate  applica¬ 
tion  this  loading  model  would  be  applied  to  penetrator  elements 
in  simplified  terradynamics  computer  codes. 

A  number  of  calculations  were  performed  using  various  pene¬ 
trator,  and  penetrator  element,  geometries  and  velocities  in 
finite  difference  wave  propagation  computer  codes  (hydrocodes) . 
Loads,  in  terms  of  applied  stresses,  were  saved  from  these  calcu¬ 
lations.  They  were  analyzed  to  develop  equations  capable  of 
predicting  the  applied  stress  levels  as  functions  of  the  pene¬ 
trator  's  geometry  and  velocity.  It  was  determined  that  loading 
functions  could  be  developed  for  the  relatively  simple  two-dimen¬ 
sional  geometries  considered.  These  functions  are  presented  in 
this  report  along  with  recommendations  for  extending  the  analysis 
to  three-dimensional  situations  and  finite  thickness  targets. 

Drag  coefficients  were  developed  from  the  loading  functions 
for  relatively  simple  geometric  shapes.  These  drag  coefficients 
are  functions  of  the  reciprocal  of  the  velocity,  so  that  there 
exists  a  maximxim  depth  of  penetration  (i.e.,  a  depth  beyond 
which  the  penetrator  cannot  advance) .  The  drag  coefficients  are 
applied  to  a  blunt  penetrator  and  comparisons  made  between  a 
trajectory  prediction  using  them  and  the  trajectory  from  a  hydro¬ 
code  calculation. 
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SECTION  II 

PRELIMINARY  CONSIDERATIONS 


The  concrete  material  behavior  model  employed  for  all  hydro- 
code  calculations  discussed  in  this  report  is  presented  in  detail 
in  Reference  1.  The  model  basically  consists  of  a  hydrostat  and 
a  yield  surface. 

The  hydrostat  defines  a  relationship  between  pressure  and 
density  and  is  valid  for  a  good  quality  concrete  with  an  unconfined 
compressive  yield  strength  of  5,000  PSI.  The  hydrostat  loading 
curves,  up  to  60  kilobars  (kb),  are  defined  by  the  following 
equations : 


P  =  144y  if  0<y<  0.0025 

P  =  0.358  +  78.62  (y-0.0025)  if  0.0025<y<  0.1 
P  =  8.0  +  130.0  (y  -  0.1)  if  0.1<y<  0.2 
P  =  21.0  +  420.0  (y  -  0.2)  if  0.2<y<0.3 


In  these  equations,  P  is  given  in  kb  and  y  is  the  excess  compres¬ 
sion  defined  by: 


y  =  p/p^,  -1 


where  p  is  ambient  concrete  density  (2.2  g/cc)  and  p  is  the 
density°at  some  compressed  state. 

The  yield  surface  was  developed  from  test  data  for  concretes 
with  varying  unconfined  compressive  strengths  and  is  a  function  of 
pressure  and  the  unconfined  strength.  The  yield  strength  is 
specified  by  the  following  equations: 


if  -0.1  f'/3<P<fy3 
c  c 

if  f'/3<P<  30f' 
c  —  c 

if  30f'  <  P 
c 


Y  =  3  (P  +  0.1fy3)/l.l 

Y  =  P  +  2/3f^ 

Y  =  30.67f' 

c 


In  these  equations  f  is  the  unconfined  compressive  strength  of  the 
concrete.  The  yield^surface  is  a  Mohr-Coulomb  type  with  a  satura¬ 
tion  level  set  at  30.67f^  (153,000  PSI  for  a  5,000  PSI  concrete). 

This  concrete  model  has  been  successfully  employed  in  pene¬ 
tration  and  breaching  calculations  and  is  actively  used  at  the 
Air  Force  Armament  Laboratory,  the  Air  Force  Weapons  Laboratory, 
the  Army  Ballistic  Research  Laboratory  and  other  locations.  The 
model  is  used  in  HULL,  TOODY,  EPIC  and  other  codes  with  the 
following  constitutive  relations. 
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d  i  ”1 

Deviatoric  stresses,  T  are  computed  from  deviatoric 
strains,  by  using  the  shear  modulus,  G. 

=  2  G 

These  deviatoric  stresses  are  limited  by  the  yield  surface  accord¬ 
ing  to  the  following  equation  (in  two-dimensional  geometry) 

Jj,  =  <  yV3 

where  is  the  second  invariant  of  the  stress  deviator  tensor. 

If  ^2  exceeds  Y^/3,  all  deviatoric  stresses- are  reset  to  lie  on 

the  yield  surface.  When  the  deviatoric  stresses  have  been  reset 
(if  necessary)  they  are  added  to  the  pressure  to  compute  the  total 
stress  as  follows: 


^ii  d  ii  „ 
T  =  T  -P 


ij  d  ij  ■  /  ■ 

T-'  =  T  1  J 


This  formulation  assumes  that  pressure  is  positive  in  compression 
while  stress  is  negative  in  compression. 


In  addition  to  these  equations  deviatoric  strains  and  stresses 
must  satisfy  the  following  relationships: 


A  =  dilatation  =  -(e^^  + 

d  ii  ii  .  , 

e  =  e  +  A/3 

d  ij  ij  .  ,  . 

11  ,  d22  I  daa 


=  0 


P  =  -1/3  (t“  +  +  T^^) 


d  ii  ii  ,  ^ 

T  =  T  +  P 

d  ij  ij  .  /  . 

1  1  d^2  2  ^  3  3  _ 


=  0 


In  one-dimensional  plane  loading,  the  strains  satisfy  the 
following  relationships: 
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0 

g.22  _  g.33  _ 

g.12  _  g.23  _  g.13  _  Q 

where  the  direction  of  loading  is  the  "1"  direction.  This  means 
that  the  deviatoric  strains  satisfy  the  equations: 

=  0 

d  22  _  d_33  _ 

since  deviatoric  stresses  are  proportional  to  deviatoric  strains, 
the  deviatoric  stresses  must  be  related  to  each  other  in  the  same 
manner  as  the  strains. 

d  22  _  d  33  _  _i^d  1  1 

These  relations  can  be  put  into  the  yield  surface  equation  for 
plastic  flow  to  obtain 

=  yV3 

or 

^  =  -  (2/3)  Y 

where  the  minus  sign  indicates  compressive  flow.  Therefore 

d^22  =  d^33  ^  Y/3 

and  the  total  stresses  in  plastic  flow  then  must  be 

=  -  (2/3)  Y-P 
t22  =  t33  =  y/3  -  P 

From  the  concrete  yield  surface  equations,  P  is  approximately  Y, 
so  the  total  stresses  can  be  written  as 

=  -(2/3)P  -  P  =  -1.67P 
.^22  ^  ^33  ^  p/3  _  p  =  -0.67P 

The  ratio  of  stresses,  then,  in  the  loading  direction  to  stresses 
in  the  other  two  normal  directions  is  given  by 

^11  _  T-ii  _  -1.67  _  ^ 

22  33  -0.67 
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These  same  final  relationships  also  apply  to  spherical 
one-dimensional  geometry.  In  this  case 

£^2  -  £33  ^  0 

but  the  other  equations  are  identical. 

In  one-dimensional  cylindrical  geometry 

e  ^  ^  ^0 
^  0 

g.33  _  ^12  _  ^23  _  ^13  _  Q 

This  leads  to  the  equation 

£11  +  £22  =  -A 

But  the  dilatation.  A,  is  usually  quite  small  so  that  the  following 
approximation  is  valid: 

e  ^  ^  ^  =  0 

or 

^  =  -^£^2  and  ^£^  ^  =  0 

therefore 

and  =  0 

Placing  these  relations  into  the  yield  surface  equation  leads  to 

=  yV3 

=  -Y/^A3 

=  -Y/^/3  -  P  =  -1.58P 
t3  3  =  0  -  P  =  -P 

=  Y/^/3  -  P  =  -0.42P 

By  definition,  the  "1"  direction  is  the  radial  direction.  The 
"2'l  direction  is  the  hoop  direction.  And  the  "3"  direction  is  an 
axial  direction  normal  to  the  plane  containing  the  cylindrical 
cross-section.  The  relations  just  developed  then  predict  that 
the  ratio  of  axial  to  radial  stress  is  given  by 
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The  relationships  between  stress  components,  which  were 
developed  above,  will  be  very  useful  in  understanding  and  analyzing 
calculational  loading  results. 

The  loading  model  developed  in  this  study  assumes  a  very  large 
shock  impedance  mismatch  between  the  concrete  and  the  penetrator . 
The  mismatch  is  required  because  the  calculations  used  to  develop 
the  loading  model  assume  no  transmission  of  shock  waves  from  the 
concrete  into  the  penetrator.  A  steel  penetrator  would  meet  this 
criterion. 
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SECTION  III 

ELEMENT  LOADING  CALCULATIONS 


The  HULL  (REF  2)  and  TOODY  (REF  3)  hydrocodes  were  used  to 
investigate  loading  on  simple  one-  and  two-dimensional  element 
shapes  impacting  5,000  PSI  unconfined  compressive  strength  concrete. 
The  simple  shapes  are  seen  in  Figure  1.  In  the  cylindrical  and 
spherical  cases  the  radii  of  curvature  were  varied  from  5  to  20  cm. 
Velocities  were  varied  from  50  to  500  m/sec.  In  the  slab  section 
calculations,  angles  between  the  surface  normal  and  the  velocity 
vector  were  varied  from  0  to  60  degrees. 

In  all  calculations  the  element  surface  velocity  was  held  con¬ 
stant  and  the  calculation  was  run  until  a  steady  stress  level  was 
attained. 

Figure  2  shows  typical  results  from  TOODY  calculations  of 
slabs.  The  top  figure  presents  stress  in  kilobars  versus  time  in 
microseconds  on  a  plane  wall,  or  slab,  being  driven  with  a  velocity 
of  300  m/sec  applied  at  an  angle  of  30  degrees  to  the  surface 
normal.  Because  of  the  hysteretic  nature  of  the  concrete  stress 
versus  strain  relationship  and  the  very  low  amount  of  viscosity 
used  in  the  TOODY  runs,  the  stress  is  initially  very  highly 
oscillatory  and  retains  an  approximately  10-percent  oscillation 
after  20  microseconds.  The  true  solution  is  shown  as  a  smooth 
curve  drawn  through  the  center  of  the  oscillations.  The  bottom 
figure  presents  the  same  information  for  a  300-m/sec  velocity 
applied  at  an  angle  of  45  degrees  to  the  surface  normal.  There  is 
a  slight  peak  in  each  smoothed  curve  followed  by  a  steady  value 
which  eventually  begins  a  slow  decay  from  free  surface  relief. 

The  steady  value  in  the  top  curve  is  14.5  kb  and  it  is  10.8  kb 
in  the  bottom  curve.  With  the  velocity  of  300  m/sec  applied 
normally  to  the  slab  it  was  found  that  the  steady  stress  was 
17.2  kb.  It  is  easily  seen  that  the  oblique  impact  values  can  be 
closely  predicted  by  multiplying  the  normal  impact  value  by  the 
cosine  of  the  velocity  angle. 

Results  from  one-dimensional  plane  and  cylindrical  calcula¬ 
tions  are  seen  in  Figure  3.  In  the  figure,  steady  stress  in  kb 
is  plotted  versus  impact  velocity  in  m/sec.  The  cylindrical 
calculations  require  separate  curves  to  represent  the  radial  and 
axial  stress  components.  As  discussed  in  Section  II  these  com¬ 
ponents  should  be  related  by  the  equation 

T  axial  -  , 

T  radial  = 

It  is  clear  in  the  figure. that  they  are  so  related.  It  is  also 
clear  in  the  figure  that  the  achieved  stress  level  is  very  nearly 
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a  linear  function  of  impact  velocity.  If  x  is  the  stress  loading 
in  the  dominant  direction,  then  we  can  write 

T  =  pcV 

where  p  is  the  concrete's  ambient  density  (2.2  g/cc) ,  V  is  the 
impact  velocity  and  c  is  an  appropriate  concrete  sound  speed. 

The  plane  and  cylindrical  one-dimensional  cases  and  the  two- 
dimensional  slab  cases  can  all  be  closely  fit  if  c  is  set  to 
2.2x10®  cm/sec.  If  V  is  also  in  cm/sec  then  x  will  be  in  units 
of  dynes/cm^ .  Since  10®  dynes/cm^  is  1  kilobar,  the  conversion 
to  kilobars  is  very  simply  made.  The  conversion  to  PSI  can  be 
made  by  noting  that  69,000  dynes/cm^  is  1  PSI. 

Table  I  lists  the  results  of  several  one-dimensional  and  slab 
geometry  cases  compared  to  the  stress  predicted  from 

X  =  2.2  X  2.2  X  10®  XVX  cos  0 

where  6  is  the  angle  between  the  velocity  vector  and  the  normal 
to  the  concrete  surface.  Most  predictions  are  within  10  percent 
and  all  are  within  23  percent.  An  even  closer  fit  could  be 
obtained  were  the  fit  made  a  function  of  the  geometry  (plane  or 
cylindrical) .  However,  it  will  be  most  useful  to  subsequent  work 
to  have  a  fit  available  which  is  not  a  function  of  geometry  and 
a  23  percent  error  is  considered  to  be  well  within  the  error  ex¬ 
pected  when  the  loading  is  applied  to  various  concrete  mixtures  and 
strengths.  Fortunately,  also,  the  radius  of  curvature  makes 
little  difference  in  cylindrical  or  spherical  geometries.  Thus 
the  one  fit  can  be  used  to  describe  the  steady  stress  expected 
for  all  cases  of  interest  to  be  encountered  by  penetrator  surface 
elements  in  a  terradynamics  code. 

To  address  peak  stresses  which  will  occur  prior  to  the  steady 
stress  level  and  any  dropoff  from  the  steady  stress  level  requires 
three-dimensional  calculations.  Of  course,  normal  impacts  of 
circular  cross-section  penetrators  in  three  dimensions  can  be 
represented  in  two-dimensional  axisymmetric  coordinates.  Several 
circular  cross-section  penetrators  were  examined  in  an  attempt  to 
understand  initial  and  late-time  stress  loading.  Figure  4  shows 
the  two  basic  geometries  examined — blunt  and  conically  nosed 
penetrators . 

The  blunt-nosed  penetrators  were  run  in  HULL  at  100,  200, 

300  and  500  m/sec.  Figure  5  presents  plots  of  stress  levels  on 
the  penetrator  noses  as  functions  of  time  for  the  specified 
velocity  levels.  All  penetrators  were  1  cm  in  radius  and  were  not 
allowed  to  slow  down  due  to  the  applied  stress.  Stress  values  in 
the  plot  were  determined  by  summing  the  force  on  all  nose  zones  and 
dividing  by  the  penetrator  ^s  cross-sectional  area  (tt)  . 
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Figure  3.  Steady  Stress  Values  from  One- 
Dimensional  Calculations 
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TABLE  I 

TYPICAL  ELEMENT  LOADING  CALCULATIONS 


CALCULATIONAL 

GEOMETRY 

1 

VELOCITY-M/SEC 

CALCULATED 

STRESS-KE 

PREDICTED 

STRESS-KB 

PERCENT 

ERROR 

PLANE 

50 

2.7 

2.4 

-10.4 

PLANE 

100 

5.2 

4.8 

-  6.9 

PLANE 

200 

10.7 

9.7 

-  9.5 

PLANE 

300 

17.2 

14.5 

-15.6 

PLANE 

500 

31.5 

24.2 

-23 

PLANE  30  degrees 

300 

14.5 

12.6 

-13.3 

PLANE  45  degrees 

300 

10.8 

10.3 

-  4.9 

CYLINDRICAL, 

R=10CM 

5C 

2.3 

2.4 

5.2 

CYLINDRICAL, 

R-IOCM 

IOC 

4.8 

4.8 

0 

CYLINDRICAL, 

R=10CM 

20C 

9.3 

9.7 

4.1 

CYLINDRICAL, 

R=10CM 

30C 

13.7 

14.5 

6.0 

CYLINDRICAL, 

R-IOCM 

500 

25.5 

24.2 

-5.1 

CYLINDRICAL , 

R=5CM 

30C 

13.3 

14.5 

8.9 

CYLINDRICAL, 

R-5CM 

50C 

25.4 

24.2 

-4.7 

SPHERICAL,  R= 

aocM 

30C 

12.9 

14.5 

-12.4 

12 
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Figure  4.  Circular  Cross-Section  Penetrator 
Geometries 
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To  understand  this  loading  it  must  be  examined  in  finer  de¬ 
tail.  Figure  6  shows  the  radial,  axial  and  hoop  stresses  as  well 
as  the  mean  stress  (-pressure)  for  time  history  station  1,  for  a 
blunt  projectile  impacting  concrete  at  300  m/sec.  Time  history 
station  1  is  on  the  projectile's  axis  of  symmetry.  Only  the  axial 
stress  component  contributes  to  slowing  down  the  penetrator.  This 
axial  component  is  initially  very  high  and  follows  the  cylindrical 
flow  equations  where  the  cylinder's  radial  direction  is  in  the 
direction  of  motion  of  the  penetrator.  This  flow  rapidly  changes 
to  a  cylindrical  flow  pattern  in  which  the  cylinder's  radius  is 
normal  to  the  direction  of  motion  of  the  penetrator.  Thus  the 
axial  stress  changes  from  the  maximum  stress  component  to  0.6 
times  the  maximum  stress  component  as  the  dominant  direction  of 
flow  changes  in  the  concrete.  And  the  maximum  stress  itself 
changes  from  a  peak  value  determined  by  setting  c  to  2.9  x  10®  ' 
cm/sec . 


With  this  background  the  axial  stress  seen  in  Figure  5  can 
now  be  understood.  The  peak  values  calculated  from 


T 


peak 


2.2  X  2.9  X  10®  xV 


decay  to  steady  values  which  are  approximately  given  by 


T 


steady 


0.6  X  2.2  X  2.2  X  10®  x  V 


where  the  0.6  factor  takes  into  account  the  fact  that  the  axial 
stress  is  not  the  dominant  direction.  At  the  time  the  steady 
stress  level  occurs,  the  concrete  is  predominantly  being  driven 
.radially  away  from  the  penetrator.  This  steady  stress  value  be¬ 
gins  ultimately  to  decay  as  relief  waves  from  the  entrance  free 
surface  begin  to  overcome  the  axial  loading. 


The  effect  of  a  conical  nose  on  peak  and  steady  stress 
levels  can  be  seen  in  Figure  7.  In  this  figure  axial  stress  is 
plotted  versus  time  for  three  penetrator  shapes  impacting  at 
300  m/sec.  The  shapes  are  the  blunt  nose,  a  45-degree  half  angle 
cone  and  a  60-degree  half  angle  cone.  Stress  is  seen  to  build 
up  more  slowly  for  the  conical  penetrators  since  the  force  on  the 
penetrator  is  divided  by  the  maximum  cross-sectional  area  to 
obtain  stress,  even  when  the  entire  area  is  not  in  contact  with 
the  concrete.  Table  II  provides  a  comparison  of  calculated  peak 
and  steady  axial  stress  values  for  several  blunt  and  conical-nosed 
penetrators.  In  this  table  the  axial  stress  for  cones  is  computed 
by  considering  the  area  actually  in  contact  with  the  concrete. 

The  equations  previously  discussed  for  peak  and  steady  loading 
provide  very  accurate  predictions  for  all  of  the  cases  studied 
when  it  is  noted  that  the  penetrator 's  velocity,  V,  must  be 
inultiplied  by  sin^G,  where  0  is  the  cone's  half  angle.  This  factor 
takes  into  account  the  fact  that  VsinG  is  the  velocity  applied  to 
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Figure  6.  Station  1  Stresses  for  a  Blunt 
Projectile  at  300  n./sec 
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Figure  7. 


Axial  Stress  at  300  m/sec  for  Various 
Penetrator  Noses 


the  concrete  and  Vsin^9  is  the  component  of  that  velocity  in  the 
axial  direction. 

The  time  required  to  decay  from  peak  to  steady  stress  levels 
can  be  expected  to  be  a  function  of  the  penetrator's  diameter. 

It  does  not  appear  to  be  a  very  strong  function  of  the  penetrator's 
velocity.  Figure  8  presents  axial  stress  versus  time  for  blunt 
penetrators  with  radii  of  1  and  2  cm.  The  time  required  to  decay 
to  steady  values  appears  to  be  a  linear  function  of  the  size  of  the 
penetrator.  Based  on  these  two  calculations  and  comparisons  with 
the  other  1-cm  radius  penetrator  loads ,  it  appears  that  the  time 
to  decay  from  peak  to  steady  stress  levels  can  be  approximately 
expressed  as 


T  =  14  xlO  ^  X  D(cm)  sec. 

Where  D  is  the  penetrator's  diam.eter,  it  also  appears  that  a 
straight  line  adequately  describes  this  drop  in  stress  level  with 
time. 


Eventually,  rarefactions  from  the  entrance  free  surface 
reduce  the  stress  loading  on  the  penetrator.  All  compressive 
waves  generated  in  the'  concrete  eventually  generate  relief  waves 
at  the  free  surface.  Geometrical  divergence  of  the  compressive 
waves  and  subsequent  divergence  of  the  relief  waves  generated  by 
them  insure  that  the  rarefactions  cannot  individually  reduce 
loading  by  any  significant  amount.  However,  after  a  relatively 
long  period  of  time  these  relief  waves  are  capable  of  almost 
totally  reducing  all  compressive  loading  on  the  penetrator. 

It  is  reasonable  to  suppose  that  the  time  required  for 
entrance  surface  relief  waves  to  significantly  affect  loading  will 
be  a  function  of  velocity.  The  penetrator  will  be  further 
from  the  surface,  at  some  given  time,  at  a  higher  velocity.  And, 
in  fact,  the  calculations  indicate  that  this  is  true.  For  blunt 
projectiles,  at  100  m/sec  velocity,  loading  is  significantly  re¬ 
duced  by  45  ysec.  This  time  rises  to  50  ysec  for  a  velocity  of 
200  m/sec,  60  ysec  at  a  velocity  of  300  m/sec  and  70  ysec  at 
500  m/sec.  These  times  can  be  expressed  in  terms  of  a  number  of 
wave  transits  from  the  free  surface  to  the  penetrator's  centerline. 
For  example,  for  the  100  m/sec  penetrator  the  initial  distance  to 
the  free  surface  is  0.5D  where  D  is  the  diameter  of  the  penetrator. 
At  45  microseconds  the  distance  to  the  free  surface  is  0.5D  plus 
the  depth  of  penetration  of  100  m/sec  x  100  cm/m  x  45  ysec  =  0.45 
cm.  The  average  distance,  therefore,  from  the  penetrator  center- 
line  to  the  free  surface  is  0.5  x  (1  cm  +  1.45  cm)  or  1.2  cm. 
Assuming  a  sound  speed  in  the  concrete  of  2.2  x  10®  cm/sec,  the 
average  distance  of  1.2  cm  represents  8.2  transits.  Similarly, 
the  200  m/sec  penetrator  begins  to  lose  steady  stress  at  7.3  tran¬ 
sits,  the  300  m/sec  penetrator  loses  stress  at  6.9  transits  and 
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PEAK  AND  STEADY  STRESSES  ON  VARIOUS  PENETRATORS 
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Figure  8.  Blunt  Projectiles  Impacting 
Concrete  at  300  m/sec 
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the  500  m/sec  significantly  loses  stress  at  5.6  transits.  The 
number  of  transits  is  a  relatively  linear  function  of  velocity 
and  is  expressed  by 

A  =  8.2  -  (V-1  X  10")  X  65  X  lO"® 

where  A  is  the  number  of  transits  and  V  is  the  velocity  in  cm/sec. 
The  time,  in  seconds,  at  which  the  steady  stress  begins  to  signi¬ 
ficantly  decay  can  then  be  expressed  as 

T  =  A  X  D/(c  -  0.5  X  A  X  V) 

where  D  is  the  diameter  of  the  projectile  in  cm,  V  is  its  velocity 
in  cm/sec  and  c  is  the  sound  speed  2.2  x  10^  cm/sec. 

For  all  of  the  blunt  projectiles  run  in  this  study  the  stress 
loss  per  unit  time  is  a  constant  0.031  kb/ysec.  Computer  time  was 
not  available  to  run  any  of  the  projectiles  to  the  point  at  which 
the  stress  no  longer  decreased.  However,  for  the  purpose  of  model 
building  it  will  be  assumed  that  loss  will  continue  until  the 
stress  level  is  reduced  to  that  which  would  occur  due  to  simply 
pushing  completely  fractured  concrete,  i.e.,  . 
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SECTION  IV 

THE  ELEMENT  LOADING  MODEL 


The  penetrator  element  loading  model  resulting  from  this 
analysis  is  seen  in  Figure  9.  The  model  is  presented  in  terms 
of  applied  stress  in  dynes/cm^  versus  time  in  seconds.  Stress 
levels  are  computed  from  the  velocity  component  normal  to  the 

element's  surface.  Stress  levels  are  connected  linearly  between 
the  initial  time  of  impact  and  times  Ti  (the  time  of  peak  stress 
clearing) ,  T2  (the  time  at  which  relief  waves  from  the  entrance 
free  surface  significantly  affect  loading)  and  T3  (the  time  at 
which  applied  stress  results  from  simply  pushing  failed  concrete) . 

The  majority  of  a  penetrator ' s  loading  will  result  from  the 
relatively  long  steady  stress  level  of  0.6  x  2.2  x  2.2  x  10®  x 
Using  this  predominant  stress  level,  drag  coefficients  for  a 
normally  impacting  projectile  can  be  computed.  If  A  is  the  pro¬ 
jectile's  cross-sectional  area,  V  its  velocity,  m  its  mass  and 
its  drag  coefficient,  then 

mxV=TxA=-pxCj^xAx  V^/2 


or 

Cp  =  2t/pv2 

substituting  the  steady  value  for  t  yields 

=  2  X  0.6  X  2.2  X  2.2  x  10®  x  V/(2.2  x  V^) 

Cj^  =  2.64  X  10®/V 

Thus  Cj^  is  an  inverse  function  of  velocity.  This  means  that  as 

the  projectile  slows  its  drag  coefficient  increases.  Figure  9  is 
a  plot  of  this  drag  coefficient  for  blunt  projectiles  versus  pro¬ 
jectile  velocity.  It  is  seen  that  Cj^  significantly  increases  as 
velocity  decreases  below  200  m/sec. 

This  is  an  interesting  function  for  Cj^  since  it  means  that 

projectiles  will  actually  eventually  stop  as  opposed  to  the  endless 
trajectory  resulting  from  a  constant  Cj^  model. 

The  drag  equation  for  blunt  projectiles  becomes 

V  =  -p  X  Cj^  X  A  X  V^/(2  X  m) 

=  2.2  X  1.2  X  2.2  X  10®  x  A  x  V/(2  X  m) 

=  2.9  X  10®  xAx  V/m 
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Figure  9.  The  Element  Loading  Model 


which  has  the  solutions 

^  ^  ^  ^(-2.9  X  10®  X  A/m)t 

''o 


and 


s  =  (V-V^)  /  (-2.9  X  10®  X  A/m) 


Where  s  is  the  depth  of  penetration/  is  the  initial  velocity 

and  V  is  the  velocity  at  time  t.  These  equations _ predict  a  finite 
depth  of  penetration  since  substitution  of  V  =  0  in  the  equation 
for  s  yields 


^f inal 


/  (2.9  X  10®  X  A/m) 


Similar  equations  can  be  developed  for  conically  nosed  penetrators 
by  simply  substituting  Vsin^0  for  V,  where  6  is  the  cone's  half 
angle . 

These  equations  were  used  to  predict  penetration  depth  at^ 
various  times  for  blunt  projectiles  run  in  HULL  with  forces  being 
allowed  to  actually  slow  the  projectiles.  Two  cases  were  run.  In 
each  case  the  projectile  weighed  500  gm  and  its  cross-section  area 
was  3.14  cm^ .  In  the  first  case  impact  occurred  at  500  m/sec 
and  HULL  predicted  a  velocity  of  393  m/sec  at  a  time  of  150  ysec 
and  a  depth  in  the  semi-infinite  concrete  of  6.57  cm.  The  simpli¬ 
fied  equations  for  velocity  and  depth  predicted  a  velocity  of 
380  m/sec  and  a  depth  of  6.25  cm  at  150  ysec.  This  is  less  than  a 
5— percent  error  in  each  quantity.  In  the  second  case,  the  projec¬ 
tile  impacted  at  200  m./sec  and  the  HULL  run  predicted  a  velocity 
of  150  m/sec  and  a  depth  of  3.31  cm  at  a  time  of  193  ysec.  The 
simple  model  predicted  a  velocity  of  141  m/sec  and  a  depth  of  3.25 
cm  at  a  time  of  193  ysec.  These  predictions  are  again  within  a 
5-percent  error. 
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SECTION  V 

RECOMMENDED  FUTURE  WORK 


Several  critical  technical  issues  should  be  resolved  prior 
to  application  of  the  suggested  loading  model  in  simplified 
terradynamics  codes.  These  issues  are  discussed  below. 

The  existence  of  the  0.6  stress  multiplication  factor  should 
be  verified  in  three-dimensional  oblique  flow  fields.  It  is  be¬ 
lieved  that  the  factor  will  be  valid  in  these  flows  since  the 
projectile  element  is  pushing  concrete  basically  radially  away 
from  itself. 

The  time  at  which  the  steady  stress  begins  to  decrease  due 
to  entrance  surface  relief  waves  should  be  investigated  in 
three-dimensional  oblique  impacts.  The  simplified  equation  used 
in  the  model  described  in  this  report  assumes  a  single  distance 
from  the  element  to  the  surface.  In  an  oblique  impact  there  are 
a  multitude  of  such  distances  associated  with  every  element. 

Effects  of  interfering  flow  fields  between  elements  must 
also  be  investigated  in  three-dimensional  impacts.  The  loading 
model  assumes  that  flow  past  neighboring  elements  has  no  influence 
on  the  element  in  question. 

The  model  assumes  a  semi-infinite  concrete  thickness.  There 
are,  of  course,  no  such  concrete  structures  in  this  world. 

Loading  in  penetration  of  a  finite  thickness  concrete  slab  will 
differ  from  semi-infinite  loading  as  soon  as  the  initial  compres¬ 
sive  wave  generated  at  impact  reflects  as  a  tensile  wave  from  the 
exit  free  surface  and  reaches  the  projectile  nose. 

Further  work  on  this  model  should  concentrate  on  the  loading 
perturbations  seen  because  of  the  concrete  slab's  finite  thickness. 
When  this  stress  decrement  is  accounted  for,  the  three-dimensional 
effects  mentioned  above  should  be  addressed. 
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ATTN:  BMDATC-M,  Mr.  P.  Boyd 
Mr.  S.  Brockway 
PO  Box  1500 
Huntsville,  AL  35807 

I  HQDA  (DAMA-ARP) 

WASH  DC  20310 

1  HQDA  (DAMA-MS) 

WASH  DC  20310 

2  Commander 

US  Army  Engineer  Waterways 
Experiment  Station 
ATTN:  Dr.  P.  Hadala 

Dr.  B.  Rohani 
PO  Box  631 

Vicksburg,  MS  39180 
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1  Commander 

US  Army  Materiel  Development 
and  Readiness  Command 
ATTN:  DRCDMD-ST 

5001  Eisenhower  Avenue 
Alexandria,  VA  22333 

10  Commander 

US  Army  armament  Research 
and  Development  Command 
ATTN:  DRDAR-TD,  Dr.  R.  Weigle 

DRDAR-LC,  Dr.  J.  Frasier 
DRDAR-SC,  Dr.  D.  Gyorog 
DRDAR-LCF,  G.  Demitrack 
DRDAR-LC A, 

G.  Randers-Pehrson 
DRDAR-SCS-M, 

R.  Kwatnoski 
DRDAR-LCU,  E.  Barrieres 
DRDAR-SCM,  Dr.  E.  Bloore 
DRDAR-TSS  (2  cys) 

Dover,  NJ  07801 

2  Director 

US  Army  ARRADCOM 
Benet  Weapons  Laboratory 
ATTN:  DRDAR-LCB-TL 

Dr.  Joseph  E.  Flaherty 
Watervliet,  NY  12189 

1  Commander 

US  Army  Armament  Materiel 
Readiness  Command 
ATTN:  DRSAR-LEP-L,  Tech  Lib 

Rock  Island,  IL  61299 

1  Commander 

US  Army  Aviation  Research 
and  Development  Command 
ATTN :  DRDAV-E 

4300  Goodfellow  Blvd. 

St.  Louis,  MO  63120 
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1  Director  1  Commander 


US  Army  Air  Mobility  Research 
and  Developmen”  Laboratory 
Ames  Research  Center 
Moffett  Field,  CA  94035 

1  Commander 

US  Army  Communications  Research 
and  Developmen"  Command 
ATTN:  DRDCO-PPA-SA 
Fort  Monmouth,  NJ  07703 

1  Commander 

US  Army  Electronics  Research 
and  Development  Command 
Technical  Support  Activity 
ATTN:  DELSD-L 
Fort  Monmouth,  NJ  07703 

1  Commander 

US  Army  Missile  Command 

ATTN:  DRSMI-R 

Redstone  Arsenal,  AL  35809 

1  Commander 

US  Army  Missile  Command 
ATTN:  DRSMI-RBL 

Redstone  Arsenal,  AL  35809 

1  Commander 

US  Army  Missile  Command 
ATTN:  DRSMI-YDL 
Redstone  Arsenal,  AL  35809 

2  Commander 

US  Army  Tank  Automotive 

Research  §  Development  Cmd 
ATTN:  DRDTA-UL 

V,  H.  Pagano 
Warren,  MI  48090 


TARAIXIOM  Tank  Automotive 
Systems  Laboratory 
ATTN:  T,  Dean 

Warren,  MI  48090 

6  Director 

US  Army  Materials  and 

Mechanics  Research  Center 
ATTN:  DRXMR-T,  Mr.  J.  Bluhm 
Mr,  J,  Mescall 
Dr,  M,  Lenoe 
R.  Shea 
F,  Quigley 
DRXMR-ATL 

Watertown ,  MA  02172 
2  Commander 

US  Army  Research  Office 
ATTN:  Dr,  E.  Saibel 
Dr,  G.  Mayer 
PO  Bok  12211 
Research  Triangle  Park 
NC  27709 

1  Director 

US  Army  TRADOC  Systems 
Analysis  Activity 
ATTN:  ATAA-SL  (Tech  Lib) 

Whie  Sands  Missile  Range 
NM  83002 

1  Office  of  Naval  Research 
Department  of  the  Navy 
ATTN:  Code  ONR  439,  N.  Perrone 
800  North  Quincy  Street 
Arlington,  VA  22217 
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3  Commander 

Naval  Air  Systems  Command 
ATTN:  AIR-604 
Washington,  DC  20360 

3  Commander 

Naval  Ordnance  Systems  Command 
Washington,  DC  20360 

2  Commander 

Naval  Air  Development 
Center,  Johnsville 
Warminster,  PA  18974 

1  Commander 

Naval  Missile  Center 
Point  Mugu,  CA  93041 

2  Commander 

Naval  Ship  Engineering  Center 
ATTN:  J.  Schell 
Tech  Lib 

Washington,  DC  20362 

1  Commander  and  Director 
David  W.  Taylor  Naval  Ship 

Research  and  Development 
Center 

Bethesda,  MD  20084 

2  Commander 

Naval  Surface  Weapons  Center 
ATTN:  Dr.  W.  G.  Soper 
Mr.  N.  Rupert 
Dahlgren,  VA  22448 

4  Commander 

Naval  Surface  Weapons  Center 
ATTN:  Dr.  S.  Fishman  (2  cys) 
Code  U-11,  J.  R.  Renzi 
R.  S.  Gross 

Silver  Spring,  MD  20910 


6  Commander 

Naval  Surface  Weapons  Center 
ATTN:  Code  R-13,  F.J.  Zerilli 
K.  Kim 
E.  T.  To ton 
M.  J.  Frankel 
Code  K-22,  F.  Stecher 
J.  M.  Etheridge 
Silver  Spring,  MD  20910 

3  Commander 

Naval  Weapons  Center 
ATTN:  Code  31804,  Mr.  M.  Smith 
Code  326,  Mr.  P.  Cordle 
Code  3261, 

Mr.  T.  Zulkoski 
China  Lake,  CA  93555 

6  Commander 

Naval  Weapons  Center 
ATTN:  Code  3181,  John  Morrow 
Code  3261, 

Mr.  C.  Johnson 
Code  3171, 

Mr.  B.  Galloway 
Code  3831, 

Mr.  M.  Backman 
Mr.  R.  E.  VanDevender, Jr . 
Dr.  O.E.R.  Heimdahl 
China  Lake,  CA  93555 

2  Director 

Naval  Research  Laboratory 
ATTN:  Dr.  C.  Sanday 
Dr.  H.  Pusey 
Washington,  DC  20375 

2  Superintendent 

Naval  Postgraduate  School 
ATTN:  Dir  of  Lib 
Dr.  R.  Ball 
Monterey,  CA  93940 
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Long  Beach  Naval  Shipyard 
ATTN:  R.  Kessler 
T.  Eto 

R.  Fernandez 
Long  Beach,  CA  90822 

HQ  USAF/SAMI 
Washington,  DC  20330 

AFIS/INOT 

Washington,  DC  20330 


OOALC/MMWMC 

Hill  AFB,  UT  84406 

HQ  TAC/DRA 

Langley  AFB,  VA  23665 
TAC/INAT 

Langley  AFB,  VA  23665 

AUL-LSE  71-249 
Maxwell  AFB,  AL  36112 


10  ADTC/DLJW  (MAJ  G.  Spitale) 
Eglin  AFB,  FL  32542 

20  ADTC/DLYV  (Mr,  j;  Collins) 
Eglin  AFB,  FL  32542 

1  AFATL/DLYV 

Eglin  AFB,  FL  32542 

1  AFATL/DLODL 

Eglin  AFB,  FL  32542 

1  AFATL/CC 

Eglin  AFB,  FL  32542 

1  AFATL/DLODR 

Eglin  AFB,  FL  32542 

1  HQ  PACAF/DOOQ 

Hickam  AFB,  HI  96853 


1 

1 

1 

1 

1 

1 


AFWAL/MLLN  (Mr.  T.  Nicholas) 
Wright-Patterson  AFB,  OH  45433 

ASD/ENESS  (S.  Johns) 
Wright-Patterson  AFB,  OH  45433 

ASD/ENFEA 

Wright-Patterson  AFB,  OH  45433 
ASD/XRP 

Wright-Patterson  AFB,  OH  45433 

HQ  USAFE/DOQ 

APO  New  York  09012 


COMIPAC/I-32 
Box  38 

Camp  H.  I.  Smith,  HI  96861 


10  Battelle  Northwest  Laboratories 
ATTN :  G .  D .  Marr 

P.  0.  Box  999 
Richland,  WA  99352 


1  HQ  PACAF/OA 

Hickam  AFB,  HI  96853 


1  Director 

Lawrence  Livermore  Laboratory 
ATTN:  Ms.  C.  Westmoreland 

P.  0.  Box  808 
Livermore,  CA  94550 
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2  Lawrence  Livermore  Laboratory 

PO  Box  808 
ATTN:  Dr.  R.  Weme 

Dr.  J.O.  Hallquist 
Livermore,  CA  94550 

6  Los  Alamos  Scientific  Laboratory 

PO  Box  1663 
ATTN:  Dr.  R.  Karpp 
Dr.  J.  Dienes 
Dr.  J.  Taylor 
Dr.  E.  Fugelso 
Dr.  D.  E.  Upham 
Dr.  R.  Keyser 
Los  Alamos,  NM  87545 

5  Sandia  Laboratories 

ATTN:  Dr.  R.  Woodfin 
Dr.  M.  Sears 
Dr.  W.  Herrmann 
Dr.  L.  Bertholf 
Dr.  A.  Chabai 
Albuquerque,  NM  87115 

1  Headquarters 

National  Aeronautics  and 
Space  Administration 
Washington,  DC  20546 

1  Jet  Propulsion  Laboratory 

4800  Oak  Grove  Drive 
ATTN:  Dr.  Ralph  Chen 
Pasadena,  CA  91102 

1  Director 

National  Aeronautics  and 
Space  Administration 
Langley  Research  Center 
Langley  Station 
Hampton,  VA  23365 

1  US  Geological  Survey 

2255  N.  Gemini  Drive 
ATTN:  Dr.  D.  Roddy 
Flagstaff,  AZ  86001 


No.  of 

Copies  Organization 

1  Aerojet  Ordnance  Company 

9236  East  Hall  Road 
Downey,  CA  90241 

1  Aeronautical  Research 

Associates  of  Princeton,  Inc. 
50  Washington  Road 
Princeton,  NJ  08540 

1  Aerospace  Corporation 

2350  E.  El  Segundo  Blvd. 

ATTN:  Mr.  L.  Rubin 
El  Segundo,  CA  90009 

1  AVCO  Systems  Division 
201  Lowell  Street 
ATTN:  Dr.  Reinecke 
Wilmington,  MA  01803 

4  Battelle  Columbus  Laboratories 

505  King  Avenue 
ATTN:  Dr.  M.  F.  Kanninen 
Dr.  G.  T.  Hahn 
Dr.  L.  E.  Hulbert 
Dr.  S.  Sampath 
Columbus,  OH  43201 

3  Boeing  Aerospace  Company 

ATTN:  Mr.  R.  G.  Blaisdell 
(M.S.  40-25) 

Dr.  N.  A.  Armstrong, 

C .  J .  Artura 
(M.S.  8C-23) 

Seattle,  WA  98124 

2  Brunswick  Corporation 
4300  Industrial  Avenue 
ATTN:  P.  S.  Chang 

R.  Grover 
Lincoln,  NE  68504 

1  Computer  Code  Consultants,  Inc. 

1680  Camino  Redondo 
ATTN:  Dr.  Wally  Johnson 
Los  Alamos,  NM  87544 
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1 


1 


2 


1 


1 


1 


1 


1 


1 


No.  of 
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Dresser  Center 
PO  Box  1407 

ATTN:  Dr.  M.S.  Chawla 
Houston,  TX  77001 

Effects  Technology,  Inc. 

5383  Hollister  Avenue 
Santa  Barbara,  CA  53111 

Firestone  Defense  Research 
and  Products 
1200  Firestone  Parkvay 
ATTN:  R.  L.  Woodall 

L.  E.  Vescelius 
Akron,  OH  44317 

FMC  Corporation 

Ordnance  Engineering  Division 

San  Jose,  CA  95114 

Ford  Aerospace  and  Coiraruni cat ions 
Corporation 
Ford  Road,  PO  Box  A 
ATTN:  L.  K.  Goodwin 

Newport  Beach,  CA  92660 

General  Atomic  Company 
PO  Box  81608 
ATTN:  R.  M.  Sullivan 
F.  H.  Ho 
S.  Kwei 

San  Diego,  CA  92138 

General  Dynamics 
PO  Box  2507 
ATTN:  J.  H.  Cuadros 

Pomona,  CA  91745 

General  Electric  Company 
Armament  Systems  Department 
Burlington,  VT  05401 

President 

General  Research  Corporation 
ATTN:  Lib 

McLean,  VA  22101 


1  Goodyear  Aerospace  Corporation 

1210  Massilon  Road 
Akron,  OH  44315 

1  H.  P.  White  Laboratory 

311^  Scarboro  Road 
Street,  MD  21154 

4  Honeywell,  Inc. 

Government  and  Aerospace 
Products  Division 
ATTN:  Mr.  J.  Blackburn 
Dr.  G.  Johnson 
Mr.  R.  Simpson 
Mr.  K.  H.  Doeringsfeld 
600  Second  Street,  NE 
Hopkins,  MN  55343 

1  Hughes  Aircraft  Corporation 

ATTr>I:  Mr.  W.  Keppel 

MS  M-5,  Bldg.  808 
Tucson,  AZ  85706 

1  International  Applied 

Physics,  Inc. 

ATTN:  Mr.  H.  F.  Swift 
7546  McEwen  Road 
Centerville,  OH  45459 

2  Kaman  Sciences  Corporation 
150C  Garden  of  the  Gods  Road 
ATTN:  Dr.  P.  Snow 

Dr.  D.  Williams 
Colcrado  Springs,  CO  80933 

1  Lockheed  Palo  Alto  Research 

Laboratory 
3251  Hanover  Street 
ATT>:  Org  5230,  Bldg.  201 
Mr.  R.  Robertson 
Pale  Alto,  CA  94394 

1  Materials  Research  Laboratory, 

Ir.c. 

1  Science  Road 
Gler.wood,  IL  60427 
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2  McDonnell-Douglas  Astro¬ 
nautics  Company 
5301  Bolsa  Avenue 
ATTN:  Dr.  L.  B.  Greszczuk 
Dr.  J.  Wall 

Huntington  Beach,  CA  92647 

1  New  Mexico  Institute  of 

Mining  and  Technology 
ATTN :  TERA  Group 
Socorro,  NM  87801 

1  Northrup  Corporation 

3901  W.  Broadway 
ATTN :  R .  L .  Ramkumar 
Hawthorne,  CA  90250 

1  Nuclear  Assurance  Corporation 

24  Executive  Park  West 
ATTN:  T.  C.  Thompson 
Atlanta,  GA  30245 

1  Pacific  Technical  Corporation 
460  Ward  Drive 

ATTN:  Dr.  F.  K.  Feldmann 
Santa  Barbara,  CA  93105 

2  Schumberger  Well  Services 
Perforating  Center 
ATTN:  J.  E.  Brooks 

J .  Brookman 
PO  Box  A 

Rosharon,  TX  77543 

1  Science  Applications,  Inc. 

101  Continental  Boulevard 
Suite  310 

El  Segundo,  CA  90245 

1  Ship  Systems,  Inc. 

11750  Sorrento  Valley  Road 
ATTN:  Dr.  G.  G.  Erickson 
San  Diego,  CA  92121 


No.  of 

Copies  Organization 

1  Systems,  Science  and  Software 
PO  Box  1620 

ATTN:  Dr.  R.  Sedgwick 
La  Jolla,  CA  92038 

2  TRW 

One  Space  Park,  Rl/2120 
ATTN:  D.  Ausherman 
M.  Bronstein 
Redondo  Beach,  CA  90277 

1  United  Technologies 

Research  Center 
438  Weir  Street 
ATTN:  P.  R.  Fitzpatrick 
Glastonbury,  CT  06033 

1  US  Steel  Corporation 

Research  Center 
125  Jamison  Lane 
Monroeville,  PA  15146 

1  VPI  5  SU 

106c  Norris  Hall 
ATTN:  Dr.  M.  P.  Kamat 
Blacksburg,  VA  24061 

2  Vought  Corporation 
PO  Box  225907 
ATTN:  Dr.  G.  Hough 

Dr.  Paul  M.  Kenner 
Dallas,  TX  75265 

1  Westinghouse,  Inc. 

PO  Box  79 

ATTN:  J.  Y.  Fan 

W.  Mifflin,  PA  15122 

1  Drexel  University 

Department  of  Mechanical  Engr. 
ATTN:  Dr.  P.  C.  Chou 
32d  and  Chestnut  Streets 
Philadelphia,  PA  19104 
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1  Forrestal  Research  Center 

Aeronautical  Eng.  Laboratory 
Princeton  University 
ATTN;  Dr.  A.  Eringen 
Princeton,  NJ  08540 

3  •  Southwest  Research  Institute 

Dept,  of  Mechanical  Sciences 
ATTN:  Dr.  U.  Lindholm 
Dr.  W.  Baker 
Dr.  R.  White 
8500  Culebra  Road 
San  Antonio,  TX  78228 

4  SRI  International 
333  Ravenswood  Avenue 
ATTN:  Dr.  L.  Seaman 

Dr.  L.  Curran 
Dr.  D.  Shockey 
Dr.  A.  L.  Florence 
Menlo  Park,  CA  94025 

1  State  University  of  New 

York  at  Stony  Brook 
Department  of  Materials 
Science  and  Engineering 
ATTN:  Dr.  H.  Herman 
Stony  Brook,  NY  11790 

2  University  of  Arizona 

Civil  Engineering  Department 
ATTN:  Dr.  D.  A.  DaDeppo 

Dr.  R.  Richard 
Tucson,  AZ  85721 

1  University  of  Arizona 

School  of  Engineering 
ATTN:  Dean  R.  Gallagher 
Tucson,  AZ  85721 

1  University  of  California 

Los  Angeles 
ATTN:  Dr.  M.  Ziv 

Los  Angeles,  CA  90024 


No.  of 

Copies  Organization 

1  University  of  California 
Department  of  Physics 
ATTN:  Dr.  Harold  Lewis 
Santa  Barbara,  CA  93106 

2  University  of  California 

College  of  Engineering 

ATTN:  Prof.  W.  Goldsmith 

Dr.  A.  G.  Evans 
Berkeley,  CA  94720 

2  University  of  Delaware 

Department  of  Mechanical 
Engineering 

ATTN;  Prof.  J.  Vinson 
Prof.  B.  Pipes 
Newark,  DE  19711 

1  University  of  Denver 

Denver  Research  Institute 
ATTN;  Mr.  R.  F.  Recht 
2390  S.  University  Blvd. 
Denver,  CO  80210 

2  University  of  Florida 

Department  of  Engineering 

Sciences 

ATTN:  Dr.  R.  L.  Sierakowski 
Dr.  L.  E.  Malvern 
Gainesville,  FL  32601 

1  University  of  Oklahoma 

School  of  Aerospace, 
Mechanical  and  Nuclear 
Engineering 
ATTN:  Dr.  C.  W.  Bert 

Norman,  OK  73019 

1  University  of  Vermont 

201  Votey  Bldg. 

ATTN :  C .  Brown 
Burlington,  VT  05405 
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Aberdeen  Proving  Ground 


Dir,  USAMSAA 
ATTN:  DRXSY-D 

DRXSY-MP,  H.  Cohen 
Cdr,  USATECOM 
•ATTN:  DRSTE-TO-F 
Dir,  USAAPG 

ATTN:  Mr.  S.  Keithely,  MTD 

Dir,  USACSL,  EA 
ATTN:  DRDAR-CLB-PA 
Bldg.  E3516 
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USER  EVALUATION  OF  REPORT 


Please  take  a  few  minutes  to  answer  the  questions  below;  tear  out 
this  sheet,  fold  as  indicated,  staple  or  tape  closed,  and  place 
in  the  mail.  Your  comments  will  provide  us  with  information  for 
improving  future  reports . 

1 ,  BRL  Report  Number _ 

2.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related 
project,  or  other  area  of  interest  for  which  report  will  be  used.) 


3.  How,  specifically,  is  the  report  being  used?  (Information 
source,  design  data  or  procedure,  management  procedure,  source  of 
ideas ,  etc . ) _ 


4.  Has  the  information  in  this  report  led  to  any  quantitative 
savings  as  far  as  man-hours/contract  dollars  saved,  operating  costs 
avoided,  efficiencies  achieved,  etc.?  If  so,  please  elaborate. 


5.  General  Comments  (Indicate  what  you  think  should  be  changed  to 
make  this  report  and  future  reports  of  this  type  more  responsive 
to  your  needs,  more  usable,  improve  readability,  etc.) _ 


6.  If  you  would  like  to  be  contacted  by  the  personnel  who  prepared 
this  report  to  raise  specific  questions  or  discuss  the  topic, 
please  fill  in  the  following  information. 

Name: 


Telephone  Number: 
Organization  Address: 


